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Tl-doped PbTe (Pbi_ x TFTe) is an anomalous superconductor with a remarkably high maximum 
T c value given its relatively low carrier concentration. Here, we present results of systematic mea- 
surements of superconducting parameters for this material, for Tl concentrations up to x — 1.4%. 
We find that it is a Type II, weak-coupled BCS superconductor in the dirty limit and discuss 
implications for the applicability of the charge Kondo model recently proposed to account for su- 
perconductivity in this system. 



I. INTRODUCTION 

Tl-doped PbTe (Pbi-^Tl^Te) is a degenerate semi- 
conductor with a small carrier concentration of ~ 10 20 
holes/cm 3 or less. However, for Tl concentrations x 
beyond a critical value x c ~ 0.3% it is observed to 
superconduct^ with a maximum T c of 1.5 K for the 
highest Tl concentrations^ comparable to more metal- 
lic systems. Furthermore, thallium is the only impurity 
known to cause superconductivity in PbTe, even though 
other impurities are able to dope to similar carrier con- 
centrations and similar densities of states. Given the 
anomalously high maximum T c value, combined with the 
unusual concentration dependence, there has been con- 
siderable discussion as to the role that the Tl impurities 
play in the superconductivity of this materiali 3 ^*^*^ 

PbTe has a rocksalt structure and has been treated 
with reasonable success using ionic models (i.e., 
Pb 2+ Te 2- )i& Thallium impurities substitute on the Pb 
site, and calculations have shown that Tl + is more sta- 
ble than Tl 3+ in the PbTe latticed This implies that Tl 
impurities will act as acceptors, and indeed Hall measure- 
ments confirm that for small doping levels the hole con- 
centration increases by one hole for every Tl impurity. 7 9 
Carrier freeze-out is not observed to the lowest temper- 
atures, indicating that the dopant atoms do not behave 
as hydrogen-like impurities due to the large static dielec- 
tric constant of the host PbTe>iS However, for concen- 
trations beyond a characteristic value the Hall number is 
observed to rise at a much slower rate with x and does 



not increase beyond ~ 10 20 cm" 



-3 1.11 



suggesting that the 



additional impurities act in a self-compensating manner. 
Significantly, within the uncertainty of these measure- 
ments, this characteristic concentration is the same as 
x c , the critical concentration required for superconduc- 
tivity. It is remarkable that as x is increased beyond x c , 
T c rises linearly over two orders of magnitude from 15 
mK for x ~ 0.3% to 1.5 K for x ~ 1.5%, while the hole 
concentration varies by less than a factor of two. 

This behavior, combined with the absence of any de- 
tectable magnetic impurities in the diamagnetic suscep- 
tibility, has been interpreted as evidence that the Tl im- 



purities are present in a mixed valence state composed 
of a mixture of both Tl + and Tl 3+ valences for x > x c . 
We recently argued that anomalies in the normal state 
resistivity of Tl-doped PbTe, that are present only for 
superconducting samples (x > x c ) and not for nonsuper- 
conducting samples (x < x c )r^ might be associated 
with a charge Kondo effect involving these degenerate Tl 
valence states. Within such a scenario, the quantum va- 
lence fluctuations associated with the Tl impurities also 
provide a possible pairing mechanism for holes in the va- 
lence band of the host PbTei 

In light of the anomalous behavior of Tl-doped PbTe 
we have investigated the superconducting properties of 
single crystal samples for a range of Tl concentrations 
up to the solubility limit of approximately 1.5%. In 
this paper, we present measurements of the heat capac- 
ity and H C 2 and the resulting estimates for coherence 
length, penetration depth, Ginzburg-Landau parameter, 
and critical fields. Our measurements show that the ma- 
terial is a Type II, weak-coupled BCS superconductor in 
the dirty limit. We discuss implications of these observa- 
tions for the charge Kondo model. 



II. SAMPLE PREPARATION AND 
EXPERIMENTAL METHODS 

Single crystals of Pbi_ x Tl x Te were grown by an un- 
seeded physical vapor transport method. Polycrystalline 
source material was synthesized by combining PbTe, Te, 
and either T^Te or elemental Tl in appropriate ratios 
and sintering at 600° C, regrinding between successive 
sintering steps. For the crystal growth, broken pieces 
of source material were sealed in an evacuated quartz 
ampoule and placed in a horizontal tube furnace held at 
750°C for 7-10 days. A small temperature gradient of ap- 
proximately l-2°C/cm allowed crystals to nucleate and 
grow at one or both of the cooler ends of the ampoule. 
Each vapor growth produced several well-formed crys- 
tals up to a few millimeters in size that could be cut and 
cleaved to prepare bars for thermodynamic and trans- 
port measurements. The thallium content was measured 
by Electron Microprobe Analysis (EMPA) using PbTe, 
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Te, and Tk>Te standards. Errors in Tl content x shown 
in subsequent figures reflect the uncertainty of the micro- 
probe method for such low dopant concentrations. The 
Tl concentration for individual samples was observed to 
be homogeneous within the uncertainty of this measure- 
ment. 

The heat capacity of single crystal samples was mea- 
sured using a thermal relaxation technique in a Quantum 
Design Physical Property Measurement System. Crys- 
tals with a mass of approximately 10-15 mg were pre- 
pared with a flat surface for good thermal contact to the 
sample platform. Measurements were made in zero ap- 
plied field and in a field H > H C 2 (typically H = 0.5-1 
T). The field was oriented at an arbitrary angle to the 
crystal axes, depending on the orientation of the flat sur- 
face. 

The upper critical field H C 2 was measured for several 
values of the Tl content x by following resistive transi- 
tions as a function of temperature for different applied 
magnetic fields. The resistivity was measured using ge- 
ometric bars cleaved from the larger as-grown crystals, 
such that the current flowed along the [100] direction 
while the magnetic field was oriented parallel to the 
equivalent [001] direction. Electrical contact was made 
using Epotek H20E silver epoxy on sputtered or evapo- 
rated gold pads and showed typical contact resistances of 
1-4 n. Resistivity measurements were made at 16 Hz and 
with current densities in the range of 25 mA/cm 2 (cor- 
responding to a current of 100 /iA for low-temperature 
measurements) to 1 A/cm 2 at higher temperatures. To 
check for heating effects, resistivity data were taken for 
different current densities and for warming and cooling 
cycles for each sample. 



III. RESULTS 

Heat capacity data for representative Tl concentra- 
tions are shown in Fig. as C p /T versus T 2 for applied 
fields that totally suppress the superconductivity. For 
all samples there is considerable curvature in the data 
even at low temperatures, presumably due to the rela- 
tively low Debye temperature Qd of PbTe. Data were 
fit to C/T = 7 + [3T 2 + 5T 4 from the base temperature 
(0.3 K) up to 1 K. From (3 = A(12tt 4 /5)i?e~ 3 , where 
R = 8.314 J/(mol K) and N = 2 for PbTe, we estimate 
€>d = 168 ±4 K for x = 0%, which is consistent with pre- 
vious reportsiiii^ Thallium substitution does not sub- 
stantially affect this value but causes a clear increase in 
the electronic contribution 7, suggesting a rapid rise in 
the density of states with x. Values of 7 obtained from 
the above fits are shown in the inset to Figure 1 as a 
function of Tl concentration x and are in broad agree- 
ment with previously published values for polycrystalline 
samples ~ 

Heat capacity data in zero field are shown in Figure 
2 for representative Tl concentrations with T c above 0.3 
K. T c values were obtained from the midpoint of the heat 
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FIG. 1: Heat capacity of Pbi-^TbTe single crystals, shown 
as Cp/T versus T 2 , for representative Tl concentrations. For 
superconducting samples, data were taken in an applied field 
H > H C 2- Inset shows electronic contribution 7 (left axis) 
and density of states at the Fermi level N(0) (right axis) as 
a function of Tl concentration x. Dashed line shows values 
calculated from known PbTe band parameters and measured 
values of the Hall number, as described in the main text. 
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FIG. 2: C p versus T in zero applied field showing the super- 
conducting anomaly for several Tl concentrations x. Dashed 
lines show the geometric construction used to obtain AC and 
T c for x — 1.4%. Inset shows AC/yT c as a function of x. 
Uncertainty in AC/jTc is derived principally from errors in 
the geometric construction used to estimate AC. 



capacity anomaly and agree well with data obtained from 
resistive transitions (Figure 3). The jump at T c , AC, can 
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FIG. 3: Superconducting transition temperatures T c as a 
function of Tl concentration x. Solid symbols are obtained 
from heat capacity measurements, and open symbols are from 
resistivity data. Line is drawn to guide the eye. 




FIG. 4: Representative resistivity data for x = 1.4%, showing 
the superconducting transition as a function of temperature 
for different magnetic fields (0, 108, 1083, 2166, 3791, and 
4875 Oe) applied parallel to the [001] direction. 



be estimated using a standard geometric construction ex- 
trapolating normal state and superconducting state be- 
haviors towards T c , as indicated by dashed lines for x 
— 1.4% in Figure 2. Resulting estimates of AC/jT c are 
shown in the inset to Figure 2 as a function of Tl concen- 
tration x. The value for the highest Tl concentration, x 
= 1.4%, is AC/-/T C = 1.55 ± 0.12, which is close to the 
weak coupling BCS result of 1.43. As x is reduced, the 
data show a small but significant systematic variation, 
tending towards a smaller value for smaller Tl concen- 
trations. The smallest value, 1.00 ± 0.20, is recorded for 
x = 0.8%, which is the lowest Tl concentration for which 
we can confidently extract AC given the base tempera- 
ture of our instrument. 

The upper critical field H&(T) was determined from 
resistivity measurements for several Tl concentrations. 
Representative data, showing the uniform suppression of 
the superconducting transition in an applied field, are 
shown in Fig. 0] for x — 1.4%. An estimate of T c was 
obtained from the midpoint of the resistive transition for 
each applied field. Resulting H C 2 curves are shown in 
Fig. El for x = 1.1% and 1.4%. Error bars indicate the 
width of the superconducting transition measured by the 
difference in temperature between 10% and 90% of the re- 
sistive transition. The upper critical field at zero temper- 
ature H C 2(T = 0) can be estimated from these data using 
the Werthamer-Helfand-Hohenberg approximation^ 



H c2 {0) = 0.69(dH c2 /dT) T T c . 



(1) 



Resulting values for x — 1.1% and 1.4% are approxi- 
mately 3900 Oe and 6000 Oe respectively (Table |j} and 
are consistent with reasonable extrapolations of the low- 
est temperature data in Fig. 5. The errors in H C 2(0) 
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FIG. 5: Temperature dependence of H C 2 for x = 1.1% and 
1.4% with the field parallel to the [001] direction. Error bars 
were determined as described in the main text. Lines are 
drawn to guide the eye. 



listed in Table [I] are estimated from the difference be- 
tween the above approximation and a parabolic fit to the 
observed data. 

Superconducting parameters such as the coherence 
length and penetration depth are dependent on the elec- 
tron mean free path I = vpr = VFfJ-m* / 'e, where vp is 
the Fermi velocity and /i is the hole mobility. From Hall 
effect measurements at 1.8 K, we find that the Hall num- 
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ber p H is - 9 x 10 19 cm- 3 for x = 1.1% and - 1 x 10 20 
for x = 1.4%iii Combining these data with measured 
values of the residual resistivity^ we find the hole mo- 
bility /i is approximately 100 cm 2 V _1 s _1 for x = 1.1% 
and 60 cm 2 V _1 s _1 for x = 1.4%. Taking into account 
the existence of both light and heavy holes at the L and 
X points in the Brillouin zone respectively^! we assume 
the elastic scattering limit holds at low temperatures such 
that II — Is = I- The Fermi level then lies in the range 
190-210 meV below the top of the valence band. Conse- 
quently, average values of vp are ~ 1.4 x 10 6 m/s for the 
L holes and 1 x 10 5 m/s for the £ holes. The resulting 
values for I are listed in Table [I] The principle contribu- 
tion to the uncertainty in this quantity arises from errors 
in the geometric factor used to calculate resistivity of 
samples. Propagation of this error is the dominant effect 
in the uncertainties of subsequent derived quantities, in- 
cluding £o and A e ff. 

The Ginzburg-Landau coherence length £(0) is calcu- 
lated from H c2 (0) by 



H c2 (0) = 



2< 2 (0)' 



(2) 



where $ = 2.0678 x 10~ 15 Tm 2 . Estimates for £(0) are 
290 A for x = 1.1% and 240 A for x = 1.4% (Table 
and should be independent of orientation for this cubic 
material. The small values of I implies that the material 
is in the dirty limit with I < £o- Therefore, the intrinsic 
coherence length £ can be extracted from the approx- 
imation £(0) ~ (^o) 1 ^ 2 , where the values are listed in 
Tabled In comparison, the BCS expression for £ is 



Co 



ahvF 



(3) 



where the BCS value of a is 0.18. Using values of £ 
derived from the dirty limit approximation, we find vf 
estimated from this formula (given in Table^J is between 
those calculated separately for the L and S holes. This 
is consistent with a mixed contribution from both carrier 
types due to the substantial scattering implied from the 
short mean free path. 

The London penetration depth for two carrier types is 
given by 



1 Hon L e 2 /i n s e 2 



\2 



rriL 



toe 



(4) 



where the superfluid densities rij, and are approxi- 
mated as the carrier densities for each carrier type, and 
rriL and to^ are the effective masses of each band. The 
corresponding values of Al are listed in Table [I] and are 
almost independent of orientation. In the dirty limit, we 
can estimate the effective penetration depth from 



A„ 



Al(V0 1/2 , 



(5) 



values of which are given in Table [I] These estimates 
are in good agreement with microwave conductivity mea- 
surements that show A(0) ~ 3 /im for x 



TABLE I: Superconducting parameters of Pbi- 
two representative Tl concentrations. 
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K 


48 ± 12 


79 ± 20 
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57 ± 14 Oe 


54 ± 13 Oe 
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Finally, we find the Ginzburg-Landau parameter using 
K = -^cff/£(0) and estimate H c and H c \ from the rela- 
tionships H c2 = V2kH c and H c \ = -^t- InK (Table 0. 



IV. DISCUSSION 

The above results indicate that Tl-doped PbTe is a 
Type II, BCS superconductor in the dirty limit, which 
is not too surprising given that the material is a doped 
semiconductor. To a large extent this observation rules 
out the possibility of more exotic scenarios for the super- 
conductivity, such as condensation of preformed pairs, at 
least for the highest Tl concentrations. Here, we discuss 
some implications for the charge Kondo model that we 
have previously proposed for this material and consider 
alternative explanations. First, we briefly reiterate the 
salient features of the charge Kondo model relevant to 
understanding these data. 

The idea of a charge Kondo effect associated with de- 
generate valence (charge) states of a valence-skipping el- 
ement was first discussed by Taraphder and Coleman, 18 
and was later re-examined in the limit of dilute impu- 
rities for the case T c ~ Tjc by Dzero and Schmalian. 7 
Weak hybridization of these degenerate impurity states 
with conduction electrons (or in the case of Tl-doped 
PbTe, with valence band holes), results in a Kondo- like 
effect with various parallels to the more common mag- 
netic case. Here, the pseudospins correspond to zero or 
double occupancy of an impurity orbital, which can be 
described in terms of a negative effective U. The degen- 
eracy of the two valence states in PbTe is not accidental 
but emerges naturally from the doping effect of the Tl 
impurities themselves^ For values of the chemical po- 
tential less than a characteristic value [i* , the impurities 
are all present as one valence (Tl + ), which act to dope 
the material. As more impurities are added, eventually a 
value of the chemical potential fx* is reached for which the 
two valence states are degenerate, and at which value the 
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chemical potential is then pinnedi The resulting charge 
Kondo effect, if present, clearly requires that hybridiza- 
tion between the impurity states and the host material is 
relatively weak. The semiconducting nature of the host 
PbTe would naturally provide an environment in which 
the local density of states at the impurity sites is rather 
small. 

Now, we discuss the origin of the enhanced electronic 
contribution to the heat capacity seen in Figure 1. The 
density of states at the Fermi energy, N(0), can be esti- 
mated from the linear term 7 in the heat capacity. The 
resulting values of N(0) are shown on the right axis of the 
inset to Figure 1. Part of the observed enhancement can 
be attributed to band filling effects, since the Hall number 
continues to rise slowly with x even for x > x c *^ How- 
ever, as has been discussed elsewhere^ the observed heat 
capacity is larger than expected from the band structure 
of PbTe alone (dashed line in the inset to Fig. 1), imply- 
ing the presence of additional states associated with the 
Tl impurities. Within the charge Kondo model, the addi- 
tional contribution would arise from the Abrikosov-Suhl 
resonance that develops at Ep for temperatures below 
Tk- If this is the case, then in principle we can esti- 
mate the concentration of Kondo impurities using the 
crude approximation jTk ~ -Rln2 per mole of impuri- 
ties. Unfortunately, uncertainty in the band parameters 
describing the £ band 21 means that it is difficult to confi- 
dently extract the magnitude of the additional contribu- 
tion to the heat capacity over and above the band-filling 
effect. Nevertheless, we can make a rough estimate to 
at least put limits on the applicability of this model. 
Using the measured Hall coefficient jJA published band 
parameters, 21 and the assumption that the band offset 
does not change with Tl doping, the band contribution 
to 7 can be estimated as shown by the dashed line in 
the inset to Fig. 1. For x = 1.4%, the observed 7 is 
0.46 mJmol _1 K -2 larger than the expected band con- 
tribution. If this enhancement is due to Kondo physics, 
then for Tk ~ 6 K (the value estimated in Ref. [I] for 
x = 0.3%), the concentration of Kondo impurities is ap- 
proximately 7x 10 18 cm -3 . In comparison, x — 1.4% cor- 
responds to a Tl concentration of 2 x 10 20 cm~ 3 . Hence, 
if a charge Kondo description is appropriate for this ma- 
terial, and if the Kondo temperature is ~ 6 K, then only 
a small fraction (~ 3%) of the Tl impurities are con- 
tributing to this effect. Within the charge Kondo model 
outlined above, this would imply that the Tl impurities 
in PbTe must be characterized by a range of fi* values, 
such that only the subset of impurities for which /1 = fj,* 
have degenerate valence states. 

There are other observations that appear to support 
this tentative conclusion. As we had previously notedji 
the magnitude of the resistivity anomaly is also less than 
what would be expected if all of the Tl impurities were 
contributing to the Kondo effect. Data for lower Tl 
concentrations, for which a reasonable fit of the low- 
temperature data can be made over an extended temper- 
ature range, indicated that approximately 1% of the Tl 



impurities contribute to the Kondo behavior i This is in 
broad agreement with the value deduced above from the 
heat capacity enhancement. In addition, the Hall num- 
ber is observed to continue to rise for x > x c ^ implying 
that the chemical potential is not pinned at one precise 
value, but rather is slowed in its progress as x increases, 
also consistent with a distribution of fi* values. 

Invoking Kondo physics of course implies a tempera- 
ture dependence to the enhancement of 7 for tempera- 
tures above Tk- Our measurements (Fig. 1) show that 
the enhancement to 7 is temperature independent for 
temperatures between 0.3 K and 1 K. However, uncer- 
tainty in this difference grows rapidly at higher tempera- 
tures due to the increasingly large phonon contribution to 
the heat capacity. As a result, we cannot unambiguously 
extract the temperature dependence of any heat capacity 
enhancement beyond the estimated Kondo temperature 
of 6 K. 

Within a BCS scenario, T c varies exponentially with 
— l/iV(0)U, where V is the pairing interaction. Figure 
6 shows ln(T c ) versus I/7 for samples with x > x c . For 
samples with T c > 0.5 K, both parameters were extracted 
from the same physical crystal. However, for samples 
with a lower critical temperature, T c was determined 
from resistivity measurements on different crystals from 
the same growth batch, introducing additional errors due 
to uncertainty in the Tl concentration. As can be seen, 
ln(T c ) scales approximately linearly with 1/7 within the 
uncertainty of the measurements. For a constant V , this 
would imply that the observed trend in T c with x (Fig. [31) 
is due to the increasing density of states (inset to Fig. 1). 
However, the situation is less clear if the charge Kondo 
picture is applicable, in which case V depends strongly 
on the Tl concentration^ and the enhancement in N(0) 
derives from Kondo physics. In the case of a supercon- 
ductor with magnetic impurities, although N(0) is en- 
hanced by this effect, the superconductivity is neverthe- 
less suppressed for T ~ Tk due to the pair-breaking effect 
associated with the rapid fluctuations in the magnetic 
moment In the case of the charge Kondo model, the 
situation is slightly more complex because the impurities 
now provide both a local pairing mechanism as well as a 
pair-breaking effect close to Tk- Consequently, the range 
of temperatures over which it is anticipated that T c will 
be suppressed is predicted to be much lower than Tk, 7 
in contrast to the case of magnetic impurities. Hence, for 
the case T c ~ Tk, the superconductivity can in principle 
benefit from the enhancement in N(0) due to the charge 
Kondo effect in a way that it cannot for magnetic impuri- 
ties. The observed trend shown in Fig. 6 may reflect this 
effect, but it is difficult to obtain quantitative estimates 
of the relative contributions to T c from the enhancement 
in N(Q) and the pairing interaction itself in this crossover 
regime of T c ~ Tici 

In the charge Kondo model, if T c is large compared to 
Tk, then the pseudospin moment is unscreened at T c , in 
which case the superconductivity is born from preformed 
pairs. In this limit, one would anticipate a much smaller 
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in this material, cannot be ignored. Nevertheless, in the 
absence of experiments directly probing the Tl valence 
(which are currently in progress), we cannot rule out less 
exciting possibilities, including the formation of a narrow 
impurity band with a relatively large density of states. In 
such a case, the pairing interaction would most likely be 
phonon mediated, though the substantial residual resis- 
tance might argue that strong Coulomb scattering also 
plays a role. The observed low-temperature resistivity 
anomaly would then presumably have its origin in some 
form of weak localization, though the temperature and 
field dependence of this feature appear to argue against 
such a scenario. 1 
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FIG. 6: Plot of ln(T c ) vs. I/7. Dashed line is a guide for the 
eye. 

anomaly in the heat capacity AC/jT c than the BCS re- 
sult of 1.43. As noted in Section ITTT1 this is clearly not 
the case for the highest Tl concentrations, consistent with 
our previous observation that T c ~ Tk for this material, 1 
However, it is difficult to understand the observed x de- 
pendence of AC/jT c within this same picture. Since T c 
decreases with decreasing x (Fig. 3), one would expect 
the superconductivity to become more BCS-like at lower 
Tl concentrations. Instead, we find that AC/jT c be- 
comes substantially smaller as x is reduced (inset to Fig- 
ure 2). Experiments are in progress to measure the heat 
capacity of samples with yet smaller Tl concentrations to 
even lower temperatures to see whether this trend con- 
tinues. 

Could the superconductivity in Tl-doped PbTe have its 
origin in more mundane physics after all? While the data 
presented here enable us to characterize this material as a 
BCS superconductor, they do not allow us to distinguish 
between different pairing mechanisms. As we have previ- 
ously arguedfi many aspects of the observed thermody- 
namic and transport properties are suggestive of charge 
Kondo physics. Moreover, the uniqueness of the Tl impu- 
rities, being the only dopant to cause superconductivity 



V. CONCLUSIONS 

In summary, we have shown that Tl-doped PbTe is a 
Type II, BCS superconductor in the dirty limit. None 
of these observations is in disagreement with the charge 
Kondo model previously described, though they do put 
some limitations on its applicability. Specifically, the rel- 
atively small enhancement of the electronic contribution 
to the heat capacity implies that if a charge Kondo de- 
scription is appropriate then only a small fraction of the 
Tl impurities can be participating in the Kondo physics. 
Within the model described in Ref . 0, this can be under- 
stood in terms of a distribution of n* values, such that 
only a subset of the Tl impurities have both valence states 
exactly degenerate for a particular value of the chemical 
potential within this range, though this has yet to be 
experimentally verified. 
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